function y = malagr( x,y )
m = length(x);
n = length(y);
s = zeros(n,n+1);
for i = 1:n
    t = [0,1];
B = 1;
    for j = 1:n
      if j~=i
          t = conv(t,[1,-x(j)]);
        B = B*(x(i)-x(j));
      end
    end
    s(i,:) = t*y(i)/B;
end
y =sum(s);
end
